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Abstract 


Organic peroxy (R0 2 ) and hydroperoxy (H0 2 ) radicals are key intermediates in the 
photochemical processes that generate ozone, secondary organic aerosol and reac- 
tive nitrogen reservoirs throughout the troposphere. In regions with ample biogenic 
5 hydrocarbons, the richness and complexity of peroxy radical chemistry presents a sig- 
nificant challenge to current-generation models, especially given the scarcity of mea- 
surements in such environments. We present peroxy radical observations acquired 
within a Ponderosa pine forest during the summer 2010 Bio-hydro-atmosphere inter- 
actions of Energy, Aerosols, Carbon, H 2 0, Organics and Nitrogen - Rocky Mountain 
10 Organic Carbon Study (BEACHON-ROCS). Total peroxy radical mixing ratios reach 
as high as 180pptv and are among the highest yet recorded. Using the comprehen- 
sive measurement suite to constrain a near-explicit 0-D box model, we investigate the 
sources, sinks and distribution of peroxy radicals below the forest canopy. The base 
chemical mechanism underestimates total peroxy radicals by as much as a factor of 
is 3. Since primary reaction partners for peroxy radicals are either measured (NO) or 
under-predicted (H0 2 and R0 2 , i.e. self-reaction), missing sources are the most likely 
explanation for this result. A close comparison of model output with observations re- 
veals at least two distinct source signatures. The first missing source, characterized 
by a sharp midday maximum and a strong dependence on solar radiation, is consis- 
20 tent with photolytic production of H0 2 . The diel profile of the second missing source 
peaks in the afternoon and suggests a process that generates R0 2 independently of 
sun-driven photochemistry, such as ozonolysis of reactive hydrocarbons. The maxi- 
mum magnitudes of these missing sources (~ 120 and 50pptvmin _1 , respectively) are 
consistent with previous observations alluding to unexpectedly intense oxidation within 
25 forests. We conclude that a similar mechanism may underlie many such observations. 
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1 Introduction 

Peroxy radicals are central components of the tropospheric radical pool. Organic per- 
oxy radicals (R0 2 ) are metastable intermediates in the oxidation of volatile organic 
compounds (VOC), while hydroperoxyl radicals (H0 2 ) are generated via photolysis of 
s carbonyl-containing VOC (e.g. formaldehyde), alkene ozonolysis and radical cycling 
reactions. R0 2 and H0 2 typically exhibit lifetimes of 1 to 1 000 s with respect to reaction 
with nitric oxide (NO) and other peroxy radicals; larger R0 2 may also undergo isomer- 
ization and/or unimolecular decomposition. Together with the hydroxyl (OH) and alkoxy 
(RO) radicals, these species comprise the RO x radical family. Rapid cycling among 
io RO x and the nitrogen oxide radicals (NO x = NO + N0 2 ) lies at the core of photochemi- 
cal mechanisms that regulate atmospheric composition and its associated impacts on 


air quality and climate. 

VOC + OH-+R0 2 (R1) 

VOC + hv — > H0 2 + products (R2) 

is VOC + 0 3 — > H0 2 + products (R3) 

R0 2 + NO -* RO + N0 2 (R4) 

H0 2 + NO -» OH + N0 2 (R5) 

RO + 0 2 -» H0 2 + oVOC (R6) 

H0 2 + H0 2 -> H 2 0 2 + 0 2 (R7) 

20 R0 2 + H0 2 -> ROOH + 0 2 (R8) 

R0 2 + R'0 2 -» products (R9) 


A number of key processes hinge upon the fate of peroxy radicals. For example, 
the conversion of NO to N0 2 via Reactions (R4) and (R5) is a critical step in tropo- 
spheric ozone formation (Thornton et al., 2002). Reactions of R0 2 with NO and N0 2 
25 may also form alkyl nitrates and peroxy nitrates, respectively, facilitating the redistribu- 
tion of pollutant precursors over regional and global scales (Moxim et al., 1996; Paulot 
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et al., 2012; Browne and Cohen, 2012). Conversely, cross-reactions of peroxy radi- 
cals (Reactions R7-R9) are responsible for termination of radical cycling in high-VOC, 
low-NO x regimes. These reactions also form hydrogen peroxide (H 2 0 2 ) and organic 
hydroperoxides (ROOH), which can induce oxidative stress in vegetation (Hewitt et al., 
5 1990). Transformations of R0 2 also generate oxidized VOC that may contribute to for- 
mation and growth of secondary organic aerosol (SOA), a major fraction of the global 
aerosol burden (Hallquist et al., 2009). SOA precursor production depends strongly on 
the degree of functionalization vs. fragmentation (Chacon-Madrid and Donahue, 2011), 
which in turn varies with the specific molecular structure of each R0 2 radical. Peroxy 
io radicals themselves can act as a major source of oxidants via Reaction (R5). Because 
this reaction shifts RO x into its more reactive form while generating N0 2 (a precursor 
for ozone and thus OH), it effectively amplifies atmospheric oxidizing capacity. In many 
environments, this reaction is the main daytime source of OH (Stone et al., 2012). 

The structure and abundance of VOC precursors shapes the detailed chemistry of 
is peroxy radicals. In the remote troposphere, methane and CO are the primary reactants. 
Here, methyl peroxy radical (CH 3 0 2 ) is the major R0 2 species and H0 2 concentrations 
typically exceed R0 2 (Cantrell et al., 2003). Near VOC sources, on the other hand, the 
R0 2 distribution can be considerably more complex and the prevailing chemistry is less 
well understood. This is particularly true in regions dominated by biogenic alkenes such 
20 as isoprene, 2-methyl-3-butene-2-ol (MBO) and monoterpenes. Numerous field studies 
have identified gaps in our understanding of photochemistry in these environments 
(Carslaw et al., 2001; DiGangi et al., 2011; Faloona et al., 2001; Griffith et al., 2012; 
Hofzumahaus et al., 2009; Lelieveld et al., 2008; Kim et al., 2013; Ren et al., 2008; 
Tan et al., 2001; Thornton et al., 2002; Whalley et al., 2011), while theoretical and 
25 laboratory work continues to reveal new reaction pathways that can significantly impact 
oxidant levels and reaction product distributions (Dillon and Crowley, 2008; Crounse 
et al., 2011, 2012; da Silva et al., 2010; Liu et al., 2012; Peeters and Muller, 2010; 
Peeters et al., 2009; Wolfe et al., 2012). Much of this work has focused on isoprene, 
which comprises one third of the global biogenic VOC emission budget (Guenther et al., 
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2012). In most instances, uncertainties in these reaction mechanisms center on the fate 
of first-generation R0 2 and on processes that control the balance of OH, H0 2 and R0 2 . 

Observations of total peroxy radicals in forested environments are relatively sparse. 
During the ROSE campaign in rural Alabama, Cantrell et al. (1992, 1993) reported 
s typical sunny mid-day R0 2 + H0 2 mixing ratios of 100 to 150pptv in an isoprene-rich 
environment under minimal anthropogenic influence (NO~ lOOpptv at noon). Maxi- 
mum values of up to 300pptv were observed on several days, and a steady-state 
model predicted an R0 2 /H0 2 ratio of ~ 1 for one prototypical day. Qi et al. (2005) 
also observed fairly high R0 2 + H0 2 mixing ratios (109 to 134pptv at midday) above 
io a Japanese mixed deciduous forest with typical NO mixing ratios less than 200 pptv, 
and it was noted that peroxy radical concentrations consistently peaked ~ 3 h after 
solar noon. At the PROPHET site in Northern Michigan, Mihele and Hastie (2003) 
measured midday R0 2 + H0 2 ranging from 20 to 60 pptv. As in Qi et al. (2005), peroxy 
radicals at PROPHET were found to peak several hours after solar noon, indicating 
is an important role for processes other than primary radical production. More recently, 
Griffith et al. (2012) reported midday mixing ratios of 20 to 50 pptv for the sum of H0 2 
and first-generation isoprene peroxy radicals at PROPHET during summer of 2008 and 
2009, consistent with earlier observations. Observations of H0 2 are more ubiquitous 
and are reviewed elsewhere (Stone et al., 2012), though some of these measurements 
20 may contain positive artifacts due to alkene-derived R0 2 (Fuchs et al., 2011). Investi- 
gations of RO x cycling are often constrained with measurements of OH and H0 2 but 
rarely include a constraint on RO z . Since R0 2 can comprise half or more of the to- 
tal peroxy radical budget, such observations are crucial for identifying and eliminating 
gaps in chemical mechanisms. 

25 We present an analysis of peroxy radical measurements obtained within a Pon- 
derosa pine forest during the 2010 Bio-hydro-atmosphere interactions of Energy, 
Aerosols, Carbon, H 2 0, Organics and Nitrogen - Rocky Mountain Organic Carbon 
Study (BEACHON-ROCS). Using the comprehensive suite of observations to constrain 
a near-explicit 0-D chemical box model, we examine the diel cycle of sources, sinks and 
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partitioning of peroxy radicals in this biogenic environment. Model underestimation of 
daytime peroxy radical concentrations leads us to consider potential missing radical 
sources. We quantify these missing processes and place them within the context of 
canonical chemistry. Analyzing the temporal profiles of missing peroxy radical mixing 
5 ratios and production rates, we identify potential novel mechanisms of radical genera- 
tion that are consistent with previous anomalous observations at this and other forests. 


2 Methods 

2.1 Field campaign 

BEACHON-ROCS took place from 1-31 August 2010 at the Manitou Forest Observa- 
io tory (39°06'00" N, 1 05°05'30" W, 2286 m above sea level). The research site is situated 
within a Ponderosa pine forest with an average canopy height of 18.5 m and negligible 
undergrowth. A leaf area index (LAI) of 3 for the tree canopy and a tree cover fraction 
of 0.38 gives a landscape average LAI of 1 .14. The closest major urban areas are Col- 
orado Springs (33 km SE) and Denver (70 km N). The site is occasionally impacted by 
is anthropogenic air masses, but prevailing winds bring relatively clean air from the south 
and southwest. Major biogenic VOC emissions at this site include 2-methyl-3-butene-2- 
ol (MBO), monoterpenes (MT) and methanol (Kaser et al., 2013a, b). Isoprene mixing 
ratios are typically less than 300 pptv, as this compound is only emitted from Ponderosa 
pine at very low rates (Kaser et al., 2013a). Further details on the site and observations 
20 can be found elsewhere (DiGangi et al., 2011, 2012; Karl et al., 2012; Kim et al., 2013; 
Kaser et al., 2013a, b). 

2.2 Peroxy radical measurements 

Two classes of peroxy radicals were measured via Peroxy Radical Chemical Ionization 
Mass Spectrometry (PeRCIMS), described in detail previously (Edwards et al., 2003; 
25 Hornbrook et al., 2011). The instrument was housed in a trailer on the forest floor with 
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the inlet protruding 0.5 m from the trailer wall at a height of 1.6 m and oriented to the 
southeast. The principle of detection involves three steps. First, ambient air is diluted 
with either N 2 or 0 2 , followed by addition of varying concentrations of NO and S0 2 . 
Reaction of H0 2 with NO generates OH via Reaction (R5), while reaction of R0 2 with 
s NO generates an RO radical that can react either with 0 2 to form H0 2 via Reaction (R6) 
or with NO to make RONO. 

RO + NO — > RONO (R10) 

Modulation of the NO/0 2 ratio controls the relative rates of Reactions (R6) and (R10) 
allowing the PeRCIMS to operate in either an H0 2 (high NO/0 2 ) or an R0 2 + H0 2 
io (low NO/0 2 ) measurement mode. The instrument switches between these two modes 
every 30 s. In the second and third steps, OH radicals react with S0 2 to generate 
sulfuric acid, which is then ionized via reaction with nitrate ions. Ions are detected via 
a custom-built quadrupole mass spectrometer. 

OH + S0 2 + M + 0 2 + 2H 2 0 H 2 S0 4 + H0 2 + H 2 0 + M (R1 1 ) 

is H 2 S0 4 + NO“ -» HSO; + HN0 3 (R12) 

Peroxy radical mixing ratios are reported as 1 min averages, with an uncertainty of 
±35 % and a detection limit of 2 pptv for data from each measurement mode. 

H0 2 measurements acquired via titration with NO can contain positive artifacts due 
to partial conversion of R0 2 , especially when the R0 2 radicals are produced via OH 
20 addition to alkenes. The primary type of R0 2 interfering in the H0 2 measurement are 
likely /?-hydroxyalkylperoxy radicals (/?R0 2 ), formed via OH addition to alkenes (Fuchs 
et al., 2011; Hornbrook et al., 2011; Whalley et al., 2013). Upon reaction with NO in 
the PeRCIMS inlet, these radicals quickly decompose to make an H0 2 radical that is 
then detected with nearly the same efficiency as ambient H0 2 . Laboratory experiments 
25 have shown that, relative to ambient H0 2 , the PeRCIMS sensitivity to isoprene-derived 
R0 2 is ~ 12% higher, while that to aromatic R0 2 is ~ 12% lower (Hornbrook et al., 
2011). Sensitivities to other /Thydroxyalkylperoxy radicals have not been tested, but 
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we assume their chemistry will be similar. For this reason, we define two quantities: 


H0 2 * = H0 2 + /3R0 2 
R0 2 * = H0 2 + R0 2 - H0 2 


( 1 ) 

( 2 ) 


5 HOj represents the “H0 2 mode” observations, which we assume to include H0 2 and 
/3R0 2 . ROj denotes the difference between total peroxy radicals and H0 2 , represent- 
ing an operationally-defined subset of the R0 2 pool. The partitioning of total R0 2 be- 
tween H0 2 and R0 2 is discussed further in Sects. 2.4 and 4. 

2.3 Other measurements 

io Additional observations used in this analysis include OH, NO, N0 2 , 0 3 , CO, PAN, 
PPN, formaldehyde, glyoxal, a suite of VOC (MBO, isoprene, cr-pinene, /3-pinene, 
limonene, camphene, non-speciated monoterpenes, acetone, methanol, benzene, 
toluene, methyl vinyl ketone, methacrolein, acetaldehyde, propanal, n-butanal, 1,3- 
butadiene), total OH reactivity, temperature, pressure, relative humidity and N0 2 pho- 
is tolysis frequencies. Table SI in the Supplement (SI) summarizes the measurement 
techniques and their associated uncertainties. It should be noted that some measure- 
ments were taken at various heights, and the heterogeneity of the forest canopy can 
lead to strong vertical concentration gradients for some reactive species (Wolfe et al., 
2011a). Furthermore, the openness of the MFO canopy can also leady to horizontal 
20 gradients in chemistry and composition. Most observations utilized here were taken 
at a height of 1.6-4m, within the trunk space of the canopy. These heights are also 
listed in Table SI. N0 2 photolysis frequencies were measured both above the canopy 
and at 2 m (on top of the instrument trailer, just above the OH inlet). We caution that 
horizontally-averaged in-canopy radiation is likely higher than represented by this mea- 
25 surement, as the canopy is relatively open and the V(N0 2 ) sensor was co-located with 
the relatively shaded OH inlet. Other photolysis frequencies, notably J( 0 3 ), are esti- 
mated by scaling measured V(N0 2 ) with clear-sky photolysis frequencies calculated 


from the Master Chemical Mechanism (MOM) v3.2 parameterization (Jenkin et al., 
1997; Saunders et al., 2003). 

VOC observations were acquired by four separate instruments, each with particular 
measurement heights, time resolution and speciation. A detailed comparison of these 
s measurements is presented elsewhere (Kaser et al., 2013b). Due to better temporal 
coverage, MBO, benzene, toluene and acetaldehyde observations are taken from the 
University of Innsbruck PTR-TOF-MS measurements at 25 m. Remaining VOC listed 
above are taken from the NCAR trace organic gas analyzer (TOGA) measurements 
(also at 25 m). Furthermore, mixing ratios of MBO and monoterpenes, which are emit- 
io ted by Ponderosa pine, are corrected for in-canopy gradients as described in the sup- 
plement. Vertical gradients of these compounds can be significant (up to 80% change 
in mixing ratios between 4 and 23 m), thus this correction is a necessary step towards 
accurately representing in-canopy reactivity. 

2.4 Model calculations 

is Detailed model calculations utilized the University of Washington Chemical Model 
(UWCM). This model incorporates a subset of the Master Chemical Mechanism (MOM) 
v3.2 (Jenkin et al., 1997; Saunders et al., 2003) with several updates and additional 
chemistry as detailed elsewhere (Wolfe and Thornton, 2011). Model setup is similar 
to that described by Kim et al. (2013) but with a few modifications. All input data are 
20 averaged to a 30 min diel cycle. Constraints include all observations listed above and 
constant mixing ratios of 1770ppbv for CH 4 and 550ppbv for H 2 . The MCM mech- 
anism subset includes all reactions from oxidation of MBO, isoprene, cr-pinene, /3- 
pinene, limonene, benzene, toluene, butadiene, acetaldehyde, propanal, n-butanal and 
methane. The simple MT mechanism of Wolfe and Thornton (2011) is used for cam- 
25 phene and unspeciated monoterpenes. Isomerization of first-generation isoprene per- 
oxy radicals (Peeters and Muller, 2010; Peeters et al., 2009) is also included using 
measured rate coefficients for isomerization (Crounse et al., 2011) and subsequent 
loss of hydroperoxyaldehydes (Wolfe et al., 2012); recall that isoprene is, however, 
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a relatively minor contributor to VOC chemistry at this site. To account for cloud and 
forest cover, the ratio of in-canopy observed and MCM-calculated U(N0 2 ) is used to 
generate a scaling factor, which is then applied to all MCM-parameterized photolysis 
frequencies. Some of our analysis below suggests that certain peroxy radicals may 
5 be sensitive to direct sunlight (Sect. 5.1); however, the attenuated light profile better 
represents photolytic sources of OH at the measurement location, and we use obser- 
vations of this radical to test proposed mechanisms. Emissions and deposition are not 
explicitly considered, but an additional first-order loss process with a lifetime of 24 h 
is given to all species to represent physical losses (deposition and advection/dilution) 
io and prevent buildup of long-lived products. The model is initialized at midnight and inte- 
grated for three days with observational constraints updated every 30 min. Three days 
is sufficient to reach a diurnal steady state, and results are shown from the third day. 
Specific model scenarios are described in the appropriate sections. Uncertainties in 
model mixing ratios are estimated from observational uncertainties as described in the 
is Supplement. 

Comparing model results with H0 2 requires some assumptions regarding the con- 
tribution of organic peroxy radicals to this measurement. To first order, we assume 
that the HOj measurement includes H0 2 and all first and second-generation /3- 
hydroxyalkylperoxy radicals produced from OH oxidation of MBO, isoprene, monoter- 
20 penes, MVK, MACR, butadiene, benzene and toluene; specific radicals used are listed 
in Table S2. 


3 Observations 

The majority of peroxy radical measurements were collected during the final two weeks 
of the campaign (16-30 August). Figure SI of the online supplement shows the full time 
25 series for peroxy radical mixing ratios, OH concentrations and meteorology. Meteoro- 
logical conditions were warm and moderately dry (average mid-day temperature and 
relative humidity of 24±2°C and 27 ±9%), with scattered clouds and occasional rain 
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after noon. Total peroxy radicals exhibited a regular diel cycle with daytime maxima of 
100 to 180pptv and nighttime minima of 0 to lOpptv, within the range of observations 
from several other forest sites (Cantrell et al., 1992, 1993; Qi et al.,2005). HOj tends to 
track total peroxy radicals, though the ratio of H0 2 to total peroxy radicals (not shown) 
s can vary significantly (from 0.4 to ~ 1) throughout the day. Daytime OH mixing ratios 
ranged from 3x 10 6 to 10x 10 6 moleccm -3 , while nighttime values were typically below 
the instrument detection limit of 5 x 10 5 moleccm -3 . Due to the regularity of diel cycles 
at this site and intermittent temporal overlap for many observations, our analysis will 
focus mainly on average diel behavior. All times discussed below refer to local solar 
io time. 

The mean diel cycle of peroxy radicals (Fig. 1) displays several interesting features. 
After sunrise, H0 2 and H0 2 + Ft0 2 rise synchronously with both OH and 0 3 photol- 
ysis frequency, consistent with photochemical sources. From 11:30 to 14:30, there is 
a sharp rise in HOj and H0 2 + F10 2 that exceeds the smoother diel cycle seen in OH 
is concentrations. This maximum, along with the brief spike at ~ 15:30, are consistent 
daily features and not an averaging artifact. In the afternoon, the decay of HOj is 
similar to that of OH, while H0 2 + R0 2 exhibits a shoulder that persists until ~ 16:30. 
Nighttime peroxy radicals are low but consistently above the PeRCIMS detection limit 
of 2 pptv, with mean mixing ratios from midnight to 05:00 of 4 ± 2 pptv and 6 ± 2 pptv for 
20 H0 2 and H0 2 + R0 2 , respectively. 

Figure 2 presents diel cycles for other key observations. Data are averaged over 
the full campaign (1-31 August). Though the degree of coverage varies for different 
species, most observations demonstrate consistent patterns from day to day. The sharp 
early morning rise in NO x may be due to entrainment from aloft during the breakup 
25 of the nocturnal boundary layer (Seok et al., 2013) or surface emission (Alaghmand 
et al., 2011). The onset of this feature is synchronous with a rapid rise in ozone, likely 
attributable to the former mechanism. Observed daytime NO mixing ratios of 100- 
150 pptv are typical for a rural continental site and are within the transition region 
between “high” and “low” NO x regimes, wherein both NO and other peroxy radicals 


31724 


Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper 


are expected to contribute significantly to total peroxy radical loss. MBO and monoter- 
penes, the primary emissions of Ponderosa pine, dominate the reactive VOC budget. 
Methanol is also abundant at this site (Kaser et al., 201 3b), but it is relatively inert. Oxi- 
dation products, including formaldehyde, glyoxal and PAN, build up throughout the day, 
5 peaking in the late afternoon and decaying at night due to deposition, thermal decom- 
position and other losses. OH reactivity, a measure of the total OH loss rate, maintains 
a fairly constant value of 6-7 s -1 during the day and rises to as much 1 4 s _1 at night. 


4 Model results 

Two initial model scenarios are considered. In the first simulation, all observations other 
io than peroxy radicals are used to constrain the model (“base”), while in a second sim- 
ulation OH is determined by the model (“ModOH”). Figure 3 compares model output 
with RO x observations for both scenarios, and additional model results are shown in 
the SI (Fig. S4). In the base case, total peroxy radicals are under-predicted throughout 
the day, with errors of up to a factor of 3 at midday. Un-constraining OH decreases 
is mid-day modeled peroxy radicals by more than 50 %. The diel cycle of modeled peroxy 
radicals closely tracks that of OH in both scenarios. The morning rise in modeled and 
measured peroxy radical mixing ratios differs by several hours, suggesting a role for 
processes other than OH-driven VOC oxidation in early morning production. Moreover, 
the model captures neither the strong mid-day maximum nor the afternoon shoulder. In 
20 the ModOH scenario, OH concentrations are under-predicted throughout the day. This 
missing OH is likely at least partly due to an under-prediction of H0 2 , which reduces 
the rate of OH production via Reaction (R5). 

The relationship between OH and peroxy radicals is somewhat complicated by pho- 
tochemical gradients and the role of vertical transport. Above the canopy, increased 
25 sunlight likely leads to more OH and thus more peroxy radicals, as demonstrated by 
a separate set of model runs constrained by above-canopy V(N0 2 ) data (Fig. S5). 
Some of these radicals will be transported into the canopy, sustaining photochemistry 
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amidst attenuated radiation. To our knowledge, there are no published data comparing 
within and above-canopy peroxy radical levels; however, detailed 1-D canopy modeling 
results predict relatively minor gradients in H0 2 and total R0 2 at other forests (Bryan 
et al., 2012; Makar et al., 1999; Wolfe and Thornton, 2011; Wolfe et al., 2011a). For 
s example, Wolfe and Thornton (201 1 ) calculate that peroxy radical mixing ratios change 
by less than 10% within a Ponderosa pine forest similar to that of BEACHON-ROCS. 
This may reflect the fact that decreased in-canopy production is balanced by downward 
mixing. Thus, model underestimates of peroxy radicals may be partly due to a miss- 
ing source from downward transport, but only inasmuch as production and destruction 
io rates vary between above and below-canopy environments. Moreover, the tight cou- 
pling of OH and H0 2 suggests that this phenomenon can be mostly accounted for by 
constraining OH to measurements. In this case, modeled peroxy radical mixing ratios 
are essentially independent of radiation (Fig. S5); in other words, OH-initiated chem- 
istry is the dominant peroxy radical source in the model, 
is Under-prediction of OH reactivity could also lead to disagreement between modeled 
and measured peroxy radicals. Missing OH reactivity is a common feature of many 
studies in regions dominated by biogenic VOC (Lou et al., 2010; Edwards et al., 2013). 
A detailed analysis of OH reactivity at BEACHON-ROCS shows that measured species 
can only account for 41 % of the total OH reactivity on average, with poorer agreement 
20 at night (Nakashima et al., 2013). In the current study, the model underestimates total 
measured OH reactivity by as much as 20 % during the day and by 40 to 50 % at night 
for both scenarios (Fig. S4). It is important to note that unmeasured VOC oxidation 
products comprise a substantial fraction (up to 45 %) of the modeled OH reactivity, and 
the abundance of these compounds is highly sensitive to the assumed dilution rate 
25 of 1 day -1 (Edwards et al., 2013). For example, increasing the “physical loss” lifetime 
(Sect. 2.4) from 24 to 48 h increases modeled OH reactivity by ~ 1 s _1 . Nonetheless, 
modeled and measured OH reactivity agree to within 20% during the day, suggesting 
that the base model adequately represents the overall rate of OH-driven R0 2 produc- 
tion. This implies that other processes must also influence the peroxy radical budget. 
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Overall, modeled H0 2 agrees somewhat better with observations than does total 
peroxy radicals, especially with OH constrained (Fig. 3c). Based on the above assump- 
tions, H0 2 comprises 60-70% of the modeled HOj, giving peak midday modeled H0 2 
mixing ratios of 25 and 12 pptv for the base and ModOH cases, respectively (Fig. S4). 
s Roughly 50-60% of total modeled R0 2 is produced from alkene oxidation (Fig. 4) and 
thus predicted to be detected as H0 2 . The remainder, which we refer to as R0 2 , can 
be compared with observations by taking the difference between measured H0 2 + R0 2 
and H0 2 (Fig. 3d). ROj is generally underestimated, especially in the afternoon where 
the “shoulder” from H0 2 + R0 2 appears as a distinct peak. This feature is an important 
io clue regarding the nature of missing R0 2 . 

The distribution of R0 2 radicals in the model closely follows that of VOC precur- 
sors (Fig. 4). MBO-derived R0 2 are the most abundant component during the day, 
while nighttime chemistry is dominated by monoterpenes, particularly /Tpinene and 
limonene. This trend occurs for two reasons. First, MBO emissions are dependent 
is upon both light and temperature (Harley et al., 1998), while monoterpene emissions at 
this site scale primarily with temperature (Kaser et al., 2013a), leading to different diel 
cycles in their concentrations (Fig. 2b). Secondly, MBO reacts almost exclusively with 
OH, while monoterpenes are reactive towards OH, 0 3 and the nitrate radical (N0 3 ). 
The latter oxidant typically only accumulates at night (Fry et al., 2013). Our model 
20 results indicate similar contributions from OH and N0 3 chemistry to monoterpene ox- 
idation in the first half of the night, with OH-driven loss prevailing in the early morning 
as N0 3 decays alongside its precursor, ozone. “Secondary” R0 2 arising from oxidation 
of unmeasured VOC comprise as much as 49% of modeled R0 2 , consistent with our 
earlier discussion of OH reactivity. H0 2 comprises 35-50% of the total peroxy rad- 
25 ical budget (Fig. 4). The diel cycle of H0 2 /(H0 2 + R0 2 ) is essentially the same for 
both model scenarios, demonstrating that absolute OH concentrations influence per- 
oxy radical abundances but not partitioning. The mid-day minimum in this ratio reflects 
the increased importance of Reaction (R8) as a sink of H0 2 (see also Fig. 6). 


31727 


In theory, model results can be used to estimate the contribution of H0 2 to HOj 
observations. Several potential approaches are discussed in the SI. All of these meth- 
ods inherently assume that model output faithfully represents the true peroxy radical 
distribution; however, without additional constraints on the nature of “missing” peroxy 
5 radicals, it is difficult to judge the reliability of this assumption. Thus, we elect to focus 
our analysis on the measured quantities, H0 2 and R0 2 . 


5 Analysis and discussion 

Discrepancies between modeled and observed peroxy radical mixing ratios indicate 
that the model is missing or misrepresenting sources and/or sinks of these species. 
10 Two features that stand out are the large midday maximum and the afternoon “shoul- 
der”, neither of which are captured by the model. Figure 5 quantifies this measurement- 
model mismatch. The magnitude of the model-measurement discrepancy should be 
interpreted with caution, as the combined uncertainties from model and observations 
leads to uncertainties of as much as a factor of 2 in the difference. Most of this un- 
is certainty is systematic in nature (e.g. accuracy of calibrations), thus we have more 
confidence in the diurnal pattern of these differences. Interestingly, most of the under- 
prediction in HOj arises from the midday maximum and the sharp afternoon satellite 
peak (Fig. 1), while R0 2 under-prediction persists throughout the day and includes 
most of the broad afternoon shoulder. These features allude to multiple mechanistic 
20 issues, and we examine each observation separately. 

5.1 Missing HO*: a photolytic HO 2 source? 

Understanding the nature of the midday H0 2 maximum requires that we first determine 
whether this peak is primarily due to changes in H0 2 or R0 2 . Given that (1) R0 2 does 
not display a similar feature, and (2) we expect most R0 2 to exhibit similar diel cycles 
25 (Fig. 4), additional H0 2 production is the simplest explanation. Production of specific 
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R0 2 that contribute to HOj, such as from photolysis of some yet-unidentified VOC, is 
also possible but less likely. This hypothesis implies several testable consequences, as 
detailed below. 

To gauge the magnitude of this putative H0 2 source, we first examine the rates of 
5 H0 2 production and loss calculated from the 0-D model (Fig. 6). Reaction of R0 2 with 
NO comprises ~ 60 % of the modeled H0 2 source, with smaller contributions from pho- 
tolysis and OH-reaction of oxidized VOC (mainly formaldehyde, HCHO). NO chemistry 
also dominates the loss of H0 2 , though model underestimates of peroxy radical con- 
centrations likely lead to an underestimate in the loss rates from reaction with H0 2 and 
io R0 2 . We can quantify the “missing” production rate by assuming that H0 2 is in steady 
state and calculating the loss rate of the missing H0 2 : 


^miss ^-miss [H0 2 ] m j ss /T|_|o moc | 


(3) 


Here, [H0 2 ] miss is the concentration of missing HOj (Fig. 5) and r H0 mod is the lifetime 
of H0 2 calculated from the base model scenario. At its peak, the missing production 
is rate is nearly double the total production rate from all known sources (Fig. 6, black 
line). For perspective, the maximum missing production rate of 102 pptmin _f is 8 times 
the H0 2 production rate from reaction of MBO-derived R0 2 with NO and 24 times 
the production rate from HCHO photolysis. Reducing the H0 2 lifetime to account for 
under-predicted reaction rates with H0 2 and R0 2 (estimated using total peroxy radical 
20 measurements and the modeled H0 2 /R0 2 ratio) increases P miss by less than 20%. 

Increased H0 2 will provide an additional source of OH. Thus, as a further test, we 
incorporate this extra H0 2 source directly into the ModOH model scenario and com- 
pare model-calculated OH with observations (Fig. 7). This modification significantly in- 
creases daytime OH concentrations. Modeled OH generally agrees with observations 
25 (to within combined uncertainties) in the morning and afternoon hours but is over- 
predicted in the early afternoon, concomitant with the sharp HOj maximum. The overall 
model-measurement agreement improves (slope = 0.25 vs 1.19, Fig. 7b), but the cor- 
relation degrades somewhat (r 2 = 0.81 vs. 0.75). These results suggest that such an 
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H0 2 radical source cannot be invoked without additional changes to the mechanism. 
Over-prediction of OH at the peak of extra H0 2 production likely indicates that some 
fraction of the missing HOj is actually R0 2 , counter to the simple assumption made 
above. Again, we caution that the magnitude of this missing H0 2 source is highly un- 
s certain. Furthermore, these results support the conclusion of Kim et al. (2013) that H0 2 
can be a major source of OH in the canopy environment. 

A closer analysis of the observational dataset provides further characterization of the 
missing source. Figure 8 shows an example of the relationship between peroxy radi- 
cals and solar radiation on 22 August; similar correlations were observed on many days 
io of the campaign. Cloud cover regularly reduces direct sunlight at this site, decreasing 
above-canopy V(N0 2 ) by factors of 2 to 4. Cloud effects on V(N0 2 ) measured near the 
ground are comparatively minor since most of the radiation here is a combination of 
scattered diffuse light and occasional sun flecks (i.e. direct sun). In-canopy V(N0 2 ) is 
a factor of 3 to 6 lower than above-canopy clear sky values (Fig. S5); we again cau- 
is tion that this is not necessarily representative of the “average” in-canopy environment 
(Sect. 2.3). HOj increases to ~ 80 pptv from 1 1 :30 to 1 4:30, consistent with the midday 
maximum in the diel average; however, HOj also decreases rapidly during periods of 
sustained radiation attenuation, down to levels similar to those observed before and 
after the maximum. R0 2 exhibits some correlation with above-canopy radiation (e.g. 
20 the troughs at hours 1 2:40 and 1 5:00), but generally the correlation is weaker and vari- 
ations in ROj are more independent of radiation. These trends are also borne out in 
the broader statistics of the full dataset (Fig. S7). The difference in the radiation depen- 
dence of HOj and ROj suggests that fast changes in HOj are not solely driven by the 
radiation dependence of OH. 

25 The simplest explanation for the observed behavior is production of peroxy radicals 
from gas-phase oxidation and/or photolysis of VOC. Additional production from OH 
chemistry is unlikely, as the model is constrained by measured OH and reproduces ob- 
served OH reactivity to within 20% during the day (Fig. S4). Only a handful of VOC are 
known to exclusively produce H0 2 during OH oxidation, notably formaldehyde (HCHO) 
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and glyoxal (HCOCHO), both of which are constrained by observations in our model. 
Gas-phase photolysis of an unidentified compound is another potential explanation, 
though the emission (or production) of such a molecule would need to match the unique 
profile of the HOj midday maximum (Fig. 5). None of the other 21 3 meteorological and 
5 chemical observations resemble this profile, thus we have no additional clues as to the 
nature of this source. H0 2 production from reactions of “missing” ROj with NO may 
also explain some of the missing H0 2 , but not the mid-day maximum (see Sect. 5.2 
and Fig. 11). 

Numerous investigations have inferred the presence of significant unidentified reac- 
io tive hydrocarbons in biogenic environments. Often this conclusion arises from discrep- 
ancies between measured and calculated OH reactivity (Di Carlo et al., 2004; Edwards 
et al., 2013; Lou et al., 2010; Sinha et al., 2010; Nolscher et al., 2012). It is still debated 
whether the missing reactivity is due to primary emissions or secondary oxidation prod- 
ucts, though this likely varies from site to site. While under-represented OH reactivity 
is could have a profound impact on peroxy radical chemistry elsewhere, we reiterate that 
this is not a viable explanation for missing peroxy radicals in the present study. Fast 
downward ozone fluxes (Goldstein et al., 2004; Kurpius and Goldstein, 2003) and high 
levels of oxidized VOC (Holzinger et al., 2005) have also been taken as evidence for 
unconventional in-canopy chemistry. Both of these findings originate from observations 
20 in a Ponderosa pine ecosystem similar to Manitou Experimental Forest, thus it is con- 
ceivable that similar chemistry is at play. Perhaps the most relevant line of evidence for 
the present study is the observation of unexpectedly large upward fluxes of HCHO dur- 
ing BEACHON-ROCS (DiGangi et al., 2011). Formaldehyde is a major product from the 
oxidation of nearly every VOC and is thus an excellent tracer for the overall efficiency 
25 of hydrocarbon degradation. DiGangi et al. (2011) demonstrated that a simple mass 
balance model incorporating known chemical and physical processes under-predicted 
the observed HCHO flux by a factor of 6. The investigators concluded that the miss- 
ing HCHO source could be attributed to oxidation of unidentified biogenic VOC and/or 
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direct emissions of HCHO from vegetation. In either case, closure of the HCHO flux 
budget required that the missing process correlate with solar radiation. 

It is possible that the missing sources of HCHO and HOj are related. The maximum 
“missing” HCHO flux of ~ 20pptvms _1 corresponds to an in-canopy HCHO production 
s rate of 65pptvmin _1 (for a canopy height of 18.5m), within the range of the missing 
H0 2 source (Fig. 6). Candidate precursors for both HCHO and H0 2 are methylperoxy, 
hydroxymethyl and /3-hydroxyalkoxy radicals, which decompose rapidly under normal 
atmospheric conditions: 

CH 3 0* + 0 2 — > HCHO + H0 2 (R13) 

io *CH 2 OH + 0 2 — > HCHO + H0 2 (R14) 

RR'C(0 , )CH 2 OH + 0 2 HCHO + H0 2 + RC(0)R' (R15) 

The lifetimes of these radicals are so short that the above reactions are often assumed 
to be instantaneous. In conventional chemical mechanisms, these radicals are inter- 
mediates of peroxy radical decomposition (mainly via reaction with NO). While known 
is photolytic sources for such molecules are thought to be minor, it is possible that photol- 
ysis of yet-unidentified VOC could simultaneously stimulate production of both HCHO 
and H0 2 , which in turn would accelerate HO x -driven photochemistry within and above 
the forest canopy. 

5.2 Missing RO*: evidence for unidentified VOC? 

20 In contrast to HOj, missing ROj mixing ratios exhibit a relatively smooth diel cycle 
(Fig. 5). This signal comprises most of the total missing peroxy radicals in the morning 
and afternoon, with a maximum at ~ 16:00. The shape of the diel profile is notably sim- 
ilar to that of several oxidation products, including HCHO, glyoxal and PAN (Fig. 2c). 
Many modeled oxidation products, including first-generation peroxides and organic ni- 
25 trates, also peak at this time (not shown). R0 2 sinks are dominated by reaction with 
NO (Fig. 9), and under-prediction of H0 2 in the base simulation likely results in under- 
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estimation of R0 2 loss via reaction with H0 2 . Thus, missing R0 2 reflects an issue with 
R0 2 sources. 

We can estimate the magnitude of the missing R0 2 source using an approach sim- 
ilar to that described for HOj (Eq. 3). This method requires calculation of the missing 
5 ROj lifetime; however, this value depends on the assumed structure of these peroxy 
radicals. Figure 9a illustrates this point for three representative peroxy radicals. The 
lifetime of CH 3 0 2 and MB0A02 (the primary R0 2 from MBO oxidation) ranges from 
30 to 60s throughout the day, except in the morning when NO concentrations spike 
(Fig. 2a). In contrast, the lifetime of the acetyl peroxy radical, CFI 3 CO 3 , is typically 
10 < 20 s. These differences arise mainly from NO reaction rate constants, which are 7.7, 
9.0 and 20 x 10 _ 12 cm 3 molec _1 s _1 at 298 K for Chl 30 2 , MB0A02 and CFI 3 CO 3 , re- 
spectively. Reaction with FI0 2 can also be an important sink for MB0A02 (Fig. 1 0a) and 
other large peroxy radicals. For comparison, we also show the concentration-weighted 
average R0 2 lifetime for all model species in the R0 2 group. Coincidentally, this lifetime 
15 is nearly identical to that of MB0A02 even though this radical is not included in R0 2 . 
Figure 9b compares the total production rate of modeled ROj with missing R0 2 pro- 
duction rates as calculated via Eq. (3). The magnitude of the missing production rate 
is similar to that of the “known” production rate except when RCTj, is assumed to have 
a lifetime comparable to CFI 3 CO 3 . In the latter instance, the source strength needed 
20 to explain missing R0 2 often exceeds the largest values estimated for missing H0 2 
(Fig. 6 ). For the other three cases, missing R0 2 production follows a diurnal pattern 
similar to its concentration profile (Fig. 5c) except in the morning, where the steady- 
state assumption may be invalid due to rapidly-changing NO concentrations. Recall 
that the absolute magnitude of this source is dependent on our estimate of missing 
25 ROj and thus is highly uncertain. 

The contrast between diurnal cycles of production rates for “known” and missing 
R0 2 (Fig. 9b) demonstrates that the processes driving missing ROj are not solely sun- 
driven (i.e. OH reaction with VOC). Figure 10 shows modeled tendencies for the same 
three representative R0 2 species (MBOA02, CH 3 0 2 and CH 3 C0 3 ). Overall radical 
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production tends to track with solar radiation, but some processes exhibit minor diurnal 
asymmetry. In particular, production of CH 3 0 2 and CH 3 C0 3 radicals from reactions of 
other R0 2 with NO maximizes slightly after noon due to both an increase in NO (Fig. 2a) 
and the buildup of R0 2 . Despite these features, OH is still the main driver for classical 
5 R0 2 production, and none of the 347 modeled R0 2 species exhibit a profile similar to 
that of the missing ROg. While amplification of a purely OH or light-dependent process 
may be able to explain missing ROj in the morning, such a source cannot explain the 
afternoon maximum. These results lead us to consider other non-OH R0 2 sources that 
may be under-represented in the model mechanism. 

10 One potential candidate for missing R0 2 is the acyl peroxy (AP) radical family. 
AP radicals are a special class of peroxy radical that can react with N0 2 to form 
a metastable acyl peroxy nitrate (APN): 

RC(0)0 2 + N0 2 RC(0)0 2 N0 2 (APN) (R16) 

Temperature controls the lifetime of APNs; for the conditions of BEACHON-ROCS (4 
15 to 29 °C), the lifetime of PAN ranges from 23 h to 24 min. APNs can act as a source 
or sink of AP radicals depending on equilibrium conditions and the strength of primary 
AP sources (LaFranchi et al., 2009). Our model predicts that PAN, the most abundant 
APN, is a net sink for acetyl peroxy radicals (Fig. 1 0c); however, the net rate of CH 3 C0 3 
loss via PAN formation is small compared to primary CH 3 C0 3 production and reaction 
20 with NO, indicating near-steady state conditions (Cleary et al., 2007; LaFranchi et al., 
2009). Moreover, the model predicts a number of additional APNs, with PAN (which is 
constrained by observations) comprising only 19-39% of the total budget. The next 
most abundant APN is the MCM species C4PAN5, a byproduct of MBO oxidation, 
at 17-22%. The shape of the missing ROj profile may imply an under-estimation in 
25 the source of AP radicals from decomposition of such APNs. Assuming steady state 
conditions for APNs, the AP radical concentration from this source alone is given by 


= WARN] 
*i8f[N0 2 ] 


(4) 
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where /r 18r and /r 18f are the reverse and forward rate constants for Reaction (R16). 
Based on this equation, errors in APN concentrations or reaction rates could lead to 
under-prediction of AP radicals. Errors in APN concentrations are likely not a viable 
explanation. Evidence from other investigations suggests that > 90 % of the total peroxy 
5 nitrate budget is comprised of only a handful of APNs, mostly PAN (Wooldridge et al., 
2010); thus, it is likely that modeled APN concentrations are already over-estimated. 
Errors in rate constants are more probable, as MCM rate constants for APN formation 
and loss are assumed equal to that of PAN, except for decomposition of PPN and 
MPAN, which follow IUPAC recommendations (Jenkin et al., 1997; Saunders et al., 
io 2003). Laboratory data from Kirchner et al. (1999) indicate that decomposition rates 
generally decrease with increasing size and decreasing electronegativity of the organic 
functional group. In contrast, IUPAC recommends decomposition rate constants for 
PPN and MPAN that are ~ 10% faster than that of PAN. Formation rate constants have 
not been measured for species other than PAN, though one study has suggested that 
is PPN formation may be 11 % slower than that of PAN (Sommariva et al., 2011). To 
completely explain missing R0 3 , we estimate that the equilibrium constant (A' 1 sf/^isr) 
for model APNs would need to decrease by more than a factor of 10, well beyond 
the likely uncertainty in this value. Moreover, model simulations show that most of the 
growth in the R0 2 pool from such a change is due not to AP radicals themselves but 
20 rather to the R0 2 products of the reaction of AP with NO. Thus, we conclude that AP 
radicals are not a major component of missing R0 2 . 

Other ROj generation mechanisms to consider include reaction of VOC with ozone 
or nitrate radical (N0 3 ). In the base model simulation, ozone chemistry contributes 10- 
20% to the daytime peroxy radical budget, while N0 3 chemistry is only significant at 
25 night. Ozonolysis of unidentified VOC has been invoked previously to explain anoma- 
lously high ozone fluxes (Goldstein et al., 2004; Hogg et al., 2007; Kurpius and Gold- 
stein, 2003), oxidation product concentrations (Holzinger et al., 2005) and sulfuric acid 
levels (Mauldin et al., 2012) in other forests. Decomposition of Criegee intermediates 
can simultaneously generate OH and R0 2 radicals, with measured yields ranging from 
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0.06 to near-unity (Aschmann et al., 2002; Atkinson and Arey, 2003; Shu and Atkin- 
son, 1994). In a detailed modeling study, Wolfe et al. (2011b) established an upper 
limit for R0 2 production from ozonolysis of “very reactive” VOC of 60pptvmin _1 , simi- 
lar to both our missing ROj source and the missing HCHO source inferred by DiGangi 
s et al. (2011). The latter study also determined that any missing VOC should exhibit 
a light-dependent emission profile similar to that of MBO. 

To test this hypothesis, we implement an additional set of reactions following the 
very reactive VOC mechanism described by Wolfe et al. (2011b). Specific reactions 
are listed in Table S3. Rate constants for initial oxidation of this hypothetical VOC are 
io assumed equal to those of /Tcaryophyllene, while reactions of the peroxy radical prod- 
ucts are assumed to be similar to those of the /3-pinene-derived radical BPINA02. As 
a modification to the original mechanism, we discriminate between R0 2 made by OH, 
0 3 and N0 3 chemistry, since we anticipate that OH-derived R0 2 would be detected as 
HOj in the PeRCIMS inlet. The yield of OH and R0 2 from ozonolysis is set to the up- 
is per limit of 0.1 recommended by Wolfe et al. (2011b). Very reactive VOC mixing ratios, 
shown in Fig. S8, are fixed to a diurnal cycle that scales with the observed flux of the 
sum of MBO and isoprene (Kaser et al., 2013a). The scaling factor of 0.23 is chosen to 
optimize model-measurement agreement for total peroxy radicals in the base scenario. 
We caution that inferred VRVOC mixing ratios depend directly on the assumed reaction 
20 rate constants and product yields - the product of which determines the R0 2 produc- 
tion rate. In other words, this calculation effectively constrains the VRVOC reactivity, as 
discussed in Wolfe et al. (201 1 b). We implement this mechanism for both the base and 
ModOH scenarios; results are shown in Fig. 11. Model-measurement agreement for 
all radicals improves markedly on incorporating very reactive VOC chemistry (compare 
25 to Fig. 3), though the unique diurnal patterns of H0 2 and R0 2 are not captured. With 
OH constrained to observations, model agreement with H0 2 improves due to (1) R0 2 
from reaction of very reactive VOC with OH and (2) increased H0 2 from reaction of 
new R0 2 radicals with NO. With OH determined by the model, however, HOj is over- 
predicted in the afternoon. This is mainly driven by excess OH co-produced with R0 2 
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during ozonolysis of very reactive VOC. On the other hand, the latter fraction of R02 
also improves agreement with ROj. Thus, the ozonolysis of unidentified hydrocarbons 
alone cannot provide closure of the RO x budget unless (1) OH is not produced with 
the same yield as R0 2 , contrary to canonical mechanisms, or (2) additional processes 
5 are invoked that affect the diel OH profile. Nonetheless, these results strongly support 
an intimate chemical link between missing peroxy radicals and other yet-unexplained 
observations (i.e. HCHO and ozone fluxes) within the canopy airspace. 


6 Conclusions 

Using the comprehensive suite of observations from the 2010 BEACHON-ROCS field 
io campaign, we have explored the detailed chemistry of peroxy radicals in a rural envi- 
ronment dominated by biogenic hydrocarbons. Total peroxy radical concentrations are 
among the highest yet reported, exceeding 100 pptv on every day of observations and 
reaching as high as 180 pptv. Box model calculations under-predict total peroxy radi- 
cals by as much as a factor of 3, indicative of missing sources. Though the PeRCIMS 
is instrument does not provide a strict segregation between H0 2 and R0 2 radicals, the 
data alludes to several distinct sources. High levels of HOj at mid-day, combined with 
a clear dependence on radiation and a lack of similar behavior in ROJ, suggest a miss- 
ing photolytic source of H0 2 . The magnitude of this missing source is highly uncertain, 
and its exact nature remains a mystery. R0 2 is also under-predicted, but the diel profile 
20 of missing ROj more closely resembles that of oxidation products and light-dependent 
VOC emissions. Implementing an additional chemical mechanism involving ozonolysis 
of putative “very reactive” VOC greatly improves model-measurement agreement with 
all radicals. With OH determined by the model, this same mechanism degrades agree- 
ment with measured OH and HOj unless the OH yield is assumed to be lower than the 
25 R0 2 yield, contrary to classical chemistry. While these results do not provide closure of 
the radical budget, they do imply that similar mechanisms may underlie both missing 
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peroxy radicals and other indicators of faster-than-expected chemistry at this and other 
forests. 

Failure to accurately represent peroxy radical chemistry in biogenic regimes will limit 
the reliability of model results and predictions. For example, the impact of peroxy radi- 
s cals on ozone production is a complex balance between the NO turnover rate and the 
production of radical reservoirs, such as alkyl nitrates and peroxy acyl nitrates. Modeled 
production of such compounds depends directly on the calculated R0 2 distribution, and 
changing the assumed yield of NO x reservoir species can alter ozone production on 
regional and larger scales (Farmer et al., 2011; Paulot et al., 2012). Inter-conversion 
io of reactive nitrogen also affects estimates of nitrogen deposition, which represents 
a critical link between atmospheric composition and ecosystem health (Sparks et al., 
2008). Peroxy radical chemistry is also intimately tied to atmospheric oxidizing capac- 
ity. In environments where H0 2 is a dominant precursor of OH, missing peroxy radical 
sources will lead to underestimation of the overall rate of oxidation. Understanding this 
is chemistry is particularly critical in high-VOC, low-NO x regimes where radical cycling 
is controlled by RO x + H0 2 reactions and, potentially, unimolecular decomposition of 
larger R0 2 (Archibald et al., 2010; Stavrakou et al., 2010). For all these reasons, it is 
imperative that we identify and eliminate remaining gaps in our understanding of chem- 
istry in such environments - especially if mechanistic shortcomings are as substantial 
20 as implied by this study. 

Closure of the RO x budget will require a comprehensive understanding of both pri- 
mary radical sources and the processes that control radical cycling in the continental 
boundary layer. Of critical importance is continued measurement of organic peroxy 
radicals, as these are too often a key missing constraint in modeling studies. Further 
25 speciation of these compounds would also greatly aid identification of errors in current 
chemical mechanisms. While it may not yet be technically feasible to segregate all in- 
dividual R0 2 , progress could be made by determining the distribution of different R0 2 
families (i.e. acyl peroxy, /3-hydroxyalkylperoxy, etc.). Measurements of radical reser- 
voirs and termination products, such as organic peroxides (ROOH) and alkyl nitrates 
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(R0N0 2 ), will also play a vital role in this regard. Observations of total ozone pro- 
duction (Cazorla and Brune, 2010) may also provide a useful check on observations of 
total peroxy radical concentrations while offering an integrated perspective on pertinent 
chemistry in biogenic environments. 


5 Supplementary material related to this article is available online at 
http://www.atmos-chem-phys-discuss. net/1 3/31 71 3/201 3/ 
acpd-1 3-31 71 3-201 3-supplement.pdf. 
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Fig. 1. Average diel cycles of total peroxy radicals (blue), HO 2 (red), RO!, (black), OH (green) 
and ozone photolysis frequency (filled gray area). Peroxy radicals and J( 0 3 ) are shown as 1 min 
means, while OH is displayed as a 30 min mean. J(0 3 ) is calculated by scaling measured 
J(N0 2 ) with the ratio of J(0 3 )IJ(N0 2 ) calculated from the MCM parameterization for clear- 
sky conditions. Nighttime OH values were typically below the instrument detection limit (5 x 
10 5 cm“ 3 ) and are thus set to half of this value. 
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Fig. 2. Average diel observations of (a) NO, N0 2 and ozone, (b) MBO, isoprene and total 
monoterpenes, (c) formaldehyde, glyoxal and PAN, and (d) OH reactivity. All data are 30 min 
means. Note that some species have been scaled to fit on a single axis, as denoted in the plot 
legends. 
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Fig. 3. Comparison of RO x observations with 0-D model results. Solid black lines with shaded 
gray areas represent observations and their associated uncertainties. Model scenarios include 
the base (solid blue line) and ModOH (dashed magenta line) simulations; in the latter case, 
OH concentrations are not constrained to observations. For comparison with HOj and R0 2 
in (c) and (d), the total modeled R0 2 is sub-divided into two groups as described in the text. 
Calculation of model uncertainties is described in the SI. 
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Fig. 4. Modeled distribution of organic peroxy radicals from the base simulation. Groups include 
methyl peroxy radical (red), first-generation peroxy radicals of MBO (yellow), isoprene (green) 
and monoterpenes (cyan), total acyi peroxy radicals (blue), first-generation radicals from ox- 
idation of other measured VOC (magenta), and “secondary” radicals resulting from oxidation 
of model-predicted, unmeasured VOC (gray). The ratio of H0 2 to total peroxy radicals is also 
shown (black diamonds). 
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Fig. 5. Difference between observed and modeled peroxy radical mixing ratios: (a) total peroxy 
radicals, (b) H0 2 and (c) ROj. Model values are taken from the base simulation. Shaded areas 
represent the combined uncertainty from observed and modeled mixing ratios. 
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Fig. 6. Rates of H0 2 production and loss calculated from the base model scenario. Rates are 
grouped according to the type of reaction. Production includes reaction of R0 2 with NO (gray), 
VOC photolysis (orange) and OH oxidation of VOC (green). Losses include reactions with NO 
(gray), H0 2 (blue), R0 2 (orange) and OH and ozone (green). The thick black line denotes the 
proposed missing H0 2 source, calculated as described in the text. 
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Fig. 7. Comparison of measured and modeled OH concentrations, (a) Diurnal profiles. Ob- 
servations and their uncertainties are shown as a black line with grey shading. The ModOH 
results (magenta) are the same as those shown in Fig. 3. In the H0 2 -constrained case (green), 
H0 2 is constrained to “missing” H0 2 as described in Sect. 5.1. The same uncertainties are 
assumed for both model scenarios, (b) Scatter plot of the same results for ModOH (magenta 
dots) and H0 2 -constrained (green asterisks) scenarios. Thins lines define uncertainties in mod- 
eled and measured values. Thick lines represent uncertainty-weighted major axis regression 
fits (ModOH: y = 0.24* + 0.18, r 2 = 0.81; ConH0 2 : y = 1.1 9* -0.05, r 2 = 0.75). 
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Fig. 8. Example of fast mid-day variations of solar radiation and peroxy radicals for day 234 (22 
August). H0 2 (red) and R0 2 (black) are plotted on the left axis, while N0 2 photolysis frequen- 
cies are plotted on the right axis for both above-canopy (light gray) and below-canopy (dark 
gray) measurements. All data are 1 min averages. 
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Fig. 9. (a) Modeled chemical lifetimes for several representative species: MBOA02 (solid blue 
line), methyl peroxy radical (dashed green line) and acetyl peroxy radical (dashed cyan line). 
The concentration-weighted average ROj lifetime is also shown (dash-dotted magenta line), 
(b) Missing ROj production rates calculated as described in the text. Also shown is the total 
production rate for all modeled ROjfthick black line). Model results are taken from the base 
simulation. 
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Fig. 10. Modeled chemical tendencies for three representative organic peroxy radicals: 
(a) MB0A02, (b) methyl peroxy radical and (c) acetyl peroxy radical. In (c), the loss rate at- 
tributed to N0 2 (green line) represents the net effect of PAN production and decomposition. 
Model results are taken from the base simulation. 
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Fig. 11. Comparison of RO x observations with 0-D model results using the very reactive VOC 
mechanism. Solid black lines with shaded gray areas represent observations and their asso- 
ciated uncertainties. Model simulations include the base (solid blue line) and ModOH (dashed 
magenta line) scenarios augmented with very reactive VOC chemistry; in the latter case, OH 
concentrations are not constrained to observations. For comparison with HOj and ROj in (c) 
and (d), the total modeled R0 2 is sub-divided into two groups as described in the text. Model 
uncertainties are excluded for clarity. 
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